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Abstract 

A method for computing electromagnetic properties of hadrons in lattice 
QCD is described and preliminary numerical results are presented. The elec- 
tromagnetic field is introduced dynamically, using a noncompact formulation. 
Employing enhanced electric charges, the dependence of the pseudoscalar me- 
son mass on the (anti)quark charges and masses can be accurately calculated. 
At j3 = 5.7 with Wilson action, the n + — 7r° splitting is found to be 4.9(3) 
MeV. Using the measured K° — K + splitting, we also find m u /md — .512(6). 
Systematic errors are discussed. 
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If a fundamental theory of quark masses ever emerges, it may be as important to resolve 
the theoretical uncertainty in the light quark masses as it is to accurately measure the 
top quark mass. Moreover, an accurate determination of the up quark mass might finally 
resolve the question of whether nature avoids the strong CP problem via a massless up 
quark. The particle data tables 0] give wide ranges for the up (2 < m u < 8 MeV) and 
down (5 < rrid < 15 MeV) quarks, while lowest order chiral perturbation theory |^, || 
gives m u /m,d — 0.57 ± 0.04. Numerical lattice calculations provide, in principle, a very 
precise way of studying the dependence of hadron masses on the lagrangian quark mass 
parameters B. However, the contribution to hadronic mass splittings within isomultiplets 
from electromagnetic (virtual photon) effects is comparable to the size of the up-down quark 
mass difference. Thus an accurate determination of the light quark masses requires the 
calculation of electromagnetic effects in the context of nonperturbative QCD dynamics. In 
this letter, we discuss a method for studying electromagnetic effects in the hadron spectrum. 
In addition to the SU(3) color gauge field, we introduce a U(l) electromagnetic field on the 
lattice which is also treated by Monte Carlo methods. The resulting SU(3)xU(l) gauge 
configurations are then analyzed by standard hadron propagator techniques. 

The small size of electromagnetic mass splittings makes their accurate determination 
by conventional lattice techniques difficult if the electromagnetic coupling is taken at its 
physical value. One of the main results of this paper is to demonstrate that calculations 
done at larger values of the quark electric charges (roughly 2 to 6 times physical values) 
lead to accurately measurable electromagnetic splittings in the light pseudoscalar meson 
spectrum, while still allowing perturbative extrapolation to physical values. 

The strategy of the calculation is as follows. Quark propagators are generated in the 
presence of background SU(3)xU(l) fields where the SU(3) component represents the usual 
gluonic gauge degrees of freedom, while the U(l) component incorporates an abelian photon 
field (with a noncompact gauge action) which interacts with quarks of specified electric 
charge. All calculations are performed in the quenched approximation and Coulomb gauge 
is used throughout for both components. Quark propagators are calculated for a variety of 
electric charges and light quark mass values. The gauge configurations were generated at 
j3 = 5.7 on a 12 3 x 24 lattice. 200 configurations each separated by 1000 Monte Carlo sweeps 
were used. In the results reported here, we have used four different values of charge given by 
e q =0, -0.4, +0.8, and -1.2 in units in which the electron charge is e — ^/47r/137 = .3028 .... 
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For each quark charge we calculate propagators for three light quark mass values in order to 
allow a chiral extrapolation. From the resulting 12 quark propagators, 144 quark-antiquark 
combinations can be formed. The meson propagators are then computed and masses for the 
78 independent states extracted. 

Once the full set of meson masses is computed, the analysis proceeds by a combination 
of chiral and QED perturbation theory. In pure QCD it is known that, in the range of 
masses considered here, the square of the pseudoscalar meson mass is quite accurately fit 
by a linear function of the bare quark masses Q. We have found that this linearity in the 
bare quark mass persists even in the presence of electromagnetism. For each of the charge 
combinations studied, the dependence of the squared meson mass on the bare quark mass is 
well described by lowest order chiral perturbation theory. Thus we write the pseudoscalar 
mass squared as 

m 2 P = A(e q , e f ) + m q B(e q , e 9 ) + m q B(e q , e q ) (1) 

where e q ,e q are the quark and antiquark charges, and m q ,m q are the bare quark masses, 
defined in terms of the Wilson hopping parameter by (ft -1 — K~ 1 )/2a. (Here a is the 
lattice spacing.) Because of the electromagnetic self-energy shift, the value of the critical 
hopping parameter must be determined independently for each quark charge. This is done 
by requiring that the mass of the neutral pseudoscalar meson vanish at discussed 
below. The results for the neutral pseudoscalars are shown in Figure 1. For the physical 
values of the quark charges, we expect that an expansion of the coefficients A and B in to 
first order in e 2 should be quite accurate. For the larger values of QED coupling that we use 
in our numerical investigation, the accuracy of first order perturbation theory is less clear: 
in fact, a good fit to all our data requires small but nonzero terms of order e 4 , corresponding 
to two- photon diagrams. Comparison of the order e 4 terms with those of order e 2 provides 
a quantitative check on the accuracy of QED perturbation theory. We have tried including 
all possible e 4 terms in the fit, but only retained those which significantly reduce the \ 2 P er 
degree of freedom. 

According to a theorem of Dashcn ffi], in the limit of vanishing quark mass, the value 
of to|, is proportional to the square of the total charge. Thus, we have also allowed the 
values of the critical hopping parameters for each of the quark charges to be fit parameters, 
requiring that the mass of the neutral mesons vanish in the chiral limit. Thus A takes the 
form (e q + e q ) 2 to order e 2 . (Order e 4 terms were not found unnecessary to fit the 
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(m q + m q) (GeV) 

Figure 1: The mass squared, Mp, (in GeV 2 ) for neutral pseudoscalar meson versus 
lattice bare quark masses m q + mq (in GeV) is shown for various quark charges 
e q = 0.0, -0.4,0.8 and -1.2. 
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data.) The coefficient B in (|l]) which parametrizes the slope of m 2 P may also be expanded in 
perturbation theory. Of the five possible e 4 terms in B^ 2 \e q , eg), only the e 4 , e^e^ and e^e? 
terms were found to improve the \ 2 ■ The coefficients in A and B, along with the four values 
of k c for the four quark charges, constitute a 12-parameter fit to the meson mass values. 

Before discussing the numerical results, we briefly describe the formulation of lattice 
QED which we have employed in these calculations. The gauge group in this case is abelian, 
and one has the choice of either a compact or noncompact formulation for the abelian gauge 
action. Lattice gauge invariance still requires a compact gauge-fermion coupling, but we are 
at liberty to employ a noncompact form of the pure photon action S cm - Then the theory 
is free in the absence of fermions, and is always in the nonconfining, massless phase. An 
important aspect of a noncompact formalism is the necessity for a gauge choice. We use QCD 
lattice configurations which have all been converted to Coulomb gauge for previous studies 
of heavy- light mesons. Coulomb gauge turns out to be both practically and conceptually 
convenient in the QED sector as well. 

For the electromagnetic action, we take 

^ Cm = ^ fJ.A n y — V„A„p) 2 (2) 

with e the bare electric coupling, n specifies a lattice site, V M the discrete lattice right- 
gradient in the fi direction and A n ^ takes on values between — oo and +00. Electromag- 
netic configurations were generated using (||) as a Boltzmann weight, subject to the linear 
Coulomb constraint 

ViA ni = (3) 

with V a lattice left-gradient operator. The action is Gaussian-distributed so it is a trivial 
matter to generate a completely independent set in momentum space, recovering the real 
space Coulomb-gauge configuration by Fast Fourier transform. We fixed the global gauge 
freedom remaining after the condition (|^) is imposed by setting the p = mode equal to 
zero for the transverse modes, and the p = mode to zero for the Coulomb modes on each 
time-slice. (This implies a specific treatment of finite volume effects which will be discussed 
below) . The resulting Coulomb gauge field A nfl is then promoted to a compact link variable 
^n™ = e ±tqArl ^ coupled to the quark field in order to describe a quark of electric charge ±ge. 
Quark propagators are then computed for propagation through the combined SU(3)xU(l) 
gauge field. 
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Table 1: Calculated shift of critical mass, Am c versus tadpole estimate for neutral 

pseudoscalar mesons with various quark charges, e q . All masses are in lattice units. 

e q k c Sm c ^ tadpole 

0.0 0.16923(3) 

-0.4 0.17130(2) 0.289(5) 0.251 
0.8 0.17763(3) 1.118(5) 0.942 
-1.2 0.18541(4) 2.063(6) 1.912 



Next we discuss the evaluation of critical hopping parameters for nonzero quark charge. 
The self energy shift induced by electromagnetic tadpole graphs may be computed pertur- 
batively. The one- loop tadpole graph is (for Wilson parameter r=l and at zero momentum 
in Coulomb gauge) 



<W?a/ = |r V{ — ~-\ (4) 

L4 fe4£^ 2 t sea 2 ' 



where fc M are the discrete lattice momentum components for a L 4 lattice and — sin(k^/2). 
This is entirely analogous to the well known QCD term Suiqcd §■ The mass shift is then 
given by the sum over multiple insertions at the same point, which exponentiates the one- 
loop graph. The usual strong QCD corrections at /3 = 5.7 are given in this approximation 
by an overall multiplicative factor of 1/(8k^ =0 ). Together this produces a shift of the critical 
inverse hopping parameter of 

Am c =(— ^ = — ^(l - e' 5mEM ) (5) 

The contribution from the conventional one loop radiative correction graph is found to be 
about one third the size of the tadpole. In Table 1, our numerical results for k c and the 
associated Am c is compared with the results using only the perturbative tadpole resummed 
result for the EM interactions^). 

For charge zero quarks, propagators were calculated at hopping parameter 0.161, 0.165, 
and 0.1667, corresponding to bare quark masses of 175, 83, and 53 MeV respectively. The 
gauge configurations are generated at (3 = 5.7, and we have taken the lattice spacing to be 
a -1 = 1.15 GeV as determined in Ref. After shifting by the improved perturbative values 
listed in Table 1, we select the same three hopping parameters for the nonzero charge quarks. 
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Because this shift turns out to be very close to the observed shift of K c , the quark masses for 
nonzero charge are nearly the same as those for zero charge. For all charge combinations, 
meson masses were extracted by a two-exponential fit to the pseudoscalar propagator over 
the time range t = 3 to 11. Smeared as well as local quark propagator sources were used 
to improve the accuracy of the ground state mesons masses extracted. Errors on each mass 
value are obtained by a single-elimination jackknife. The resulting data is fitted to the 
chiral/QED perturbative formula ([!]) by x 2 minimization. The fitted parameters are given 
in Table 2. Errors were obtained by performing the fit on each jackknifed subensemble. 

Aside from very small corrections of order (md — m M ) 2 , the ir + — 7r° mass splitting is 
of purely electromagnetic origin, and thus should be directly calculable by our method. 
Because we have used the quenched approximation, uu and dd mesons do not mix. The 
neutral pion mass is obtained by averaging the squared masses of the uu and dd states. (In 
full QCD the uu and dd mix in such a way that the neutral octet state remains a Goldstone 
boson of approximate chiral SU(3)xSU(3). By averaging the squared masses of uu and dd 
in the quenched calculation, we respect the chiral symmetry expected from the full theory. 
By contrast, linear averaging of the masses would give a 7r° mass squared nonanalytic in 
the quark masses). Thus, to zeroth order in e 2 , the terms proportional to quark mass Q 
cancel in the difference mz. + — m^o- This difference is then given quite accurately by the 
single term 

<»«#e 2 (6) 



rttt- 



Using the coefficients listed in Table 2, and the experimental values of the 7r°, K°, and K + 
masses, we may directly solve the resulting three equations for the up, down, and strange 
masses. The ir + — ir° splitting may then be calculated, including the very small contributions 
from the order e 2 m q terms. We obtain 



tr, 7T + - m T o = 4.9 ± 0.3MeV (7) 
compared to the experimental value of 4.6 MeV. (The electromagnetic contribution to this 



splitting is estimated 1 1 to be 4.43 ± 0.03 MeV.) Our calculation can be compared to the 
value 4.4 MeV (for Aq C d = 0.3 GeV and m s = 120 MeV) obtained by Bardeen, Bijnens 
and Gerard[[ll| using large N methods. The values obtained for the bare quark masses are 

m u = 3.86(3), m d = 7.54(5), m s = 147(1) (8) 
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The errors quoted are statistical only, and are computed by a standard jackknife procedure. 
The extremely small statistical errors reflect the accuracy of the pseudoscalar mass determi- 
nations, and should facilitate the future study of systematic errors (primarily finite volume, 
continuum extrapolation |l3| and quark loop effects), which are expected to be considerably 
larger. The relationship between lattice bare quark masses and the familiar current quark 
masses in the MS continuum regularization is perturbatively calculable |l4| . 

The presence of massless, unconfined degrees of freedom implies that the finite volume 
effects in the presence of electromagnetism may be much larger than for pure QCD. In 
fact, the corrections are expected to fall as inverse powers of the lattice size, instead of 
exponentially. We have estimated the size of the finite volume correction phenomenologically 
by considering the discussion of Bardeen, et.al|ll]], which models the low-q 2 contribution to 
the 7r + — 7T° splitting in terms of tt, p, and Al intermediate states. This gives the splitting 
as an integral, 

~ 16^2 J Q (q2 + m 2 ){q 2 + ^ d 1 W 

If the upper limit M 2 is taken to infinity, this reproduces the result of Ref . |Q , which gives 
8m n — b.lAIeV. Even better agreement with experiment is obtained by matching the low- 
q 2 behavior with the large-g 2 behavior from large N perturbative QCDjnJ. Here we only 
use the expression to estimate the finite volume correction, for which the low-g 2 expression 
above should be adequate. To estimate the finite volume effect, we cast this expression 
as a four-dimensional integral over d q and then construct the finite volume version of it 
by replacing the integrals with discrete sums (excluding the q = mode). For a 12 3 x 24 
box with a -1 = 1.15 GeV, we find that the infinite volume value of 5.1 McV is changed to 
6m n — 4.8 MeV, indicating that the result we have obtained in our lattice calculation should 
be corrected upward by about 0.3 MeV, or about 6%. In further numerical studies, we will 
be able to determine the accuracy of this estimate directly by calculations on larger box 
sizes. A study of other systematics such as finite lattice spacing effects is also in progress, 
and will be reported in a subsequent publication. 

For comparison with other results,^, |§|, || we quote the following mass ratios, which are 
independent of renormalization prescription, 

m d -m u = Q249(3) = .512(6) (10) 

m s m d 

With the errors shown, which are statistical only, these results differ significantly from 
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Table 2: Coefficients of fitting function, Eq.(|l|). Terms of order e q e^ and e 4 - in 

and e 4 in A were consistent with zero and dropped from this fit. Numerical values 

are in GeV 2 and GeV for A and B terms respectively. 

Parameter Fit 

A 0.0143(10)(e 9 + e( j) 2 

£K°) 1.594(11) 

B« 0.205(22) e 2 + 0.071(9)e 9 eg + 0.050(7)e| 

BP) 0.064(17)e^ + 0.033(6)e 3 q eq - 0.031{A)e 2 q e 2 q 

the lowest order estimate [|| which uses Dashen's theorem to estimate the electromagnetic 
contribution to the kaon splitting to zeroth order. This lowest order estimate neglects the 
quark mass dependence of the electromagnetic terms, which we have determined by our 
procedure. Specifically, the important corrections to the lowest order result come from 
terms involving the strange quark mass times the difference of up and down quark charges. 
These corrections are determined by the second and third terms in in Table 2. The 
Weinberg analysis predicts that the 4.0 MeV kaon splitting consists of 5.3 MeV from the 
up-down mass difference and -1.3 MeV from EM. In our results, the up-down mass difference 
contributes 5.9 MeV, with -1.9 MeV from EM. This goes in the direction indicated by the 
rj — > 37r decay rate pj, although our results do not deviate as much from the lowest order 
analysis as those of Ref. where the quark mass contribution to the kaon splitting is 
estimated to be 7.0 MeV. 

In the present work we have focused on the pseudoscalar meson masses. This is the 
most precise way of determining the quark masses as well as providing an important test 
of the method in the 7r + — 7r° splitting. Further calculations of electromagnetic splittings 
in the vector mesons and the baryons, as well as in heavy-light systems, are possible using 
the present method. This will provide an extensive opportunity to test the precision of the 
method and gain confidence in the results. Further study of electromagnetic properties of 
hadrons in lattice QCD, such as magnetic moments and form factors, is also anticipated. 
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